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ABSTRACT 

G54. 1+0.3 is a Crab-like pulsar wind nebula (PWN) with the highest 7-ray to X-ray 
luminosity ratio among all the nebulae driven by young rotation-powered pulsars. We 
model the spectral evolution of the PWN and find it difficult to match the observed 
multi-band data with leptons alone using reasonable model parameters. In lepton- 
hadron hybrid model instead, TeV photons come mainly from ir° decay in proton- 
proton interaction and the observed photon spectrum can be well reproduced. The 
newly discovered infrared loop and molecular cloud in or closely around the PWN can 
work as the target for the bombardment of the PWN protons. 

Key words: gamma rays: theory - ISM: individual (G54. 1+0.3) - radiation mecha- 
nisms: non-thermal 



1 INTRODUCTION 

Pulsar wind nebulae (PWNe) are thought to be an efficient 
accelerator for cosmic rays with energy above the "knee". 
Pulsar, located in the center of PWN, loses its energy by 
driving ultra-relativistic wind of electrons, positrons, and 
ions. However, it is hard to know the fraction of energy di- 
vision of different particle components. The extended 7-ray 
emission from PWN provides an exciting opportunity for 
studying the acceleration and radiation mechanism of par- 
ticles in ultra-relativistic shocks. It has been long debated 
whether the very high energy (VHE) emission from PWNe 
as well as from supernova remnants (SNRs) is leptonic or 
hadronic origin. Theoretically, it has been suggested that 
some fraction of the pulsar's spin-down energy can be con- 
verted into nuclei (Cheng et al. 1990; Arons & Tavani 1994), 
which indicates that TeV emission from PWNe may con- 
tain contribution from both leptons and hadrons. Indeed, 
nucleonic models have been used to reproduce the 7-rays 
from Crab and Vela X, respectively (Atoyan et al. 1996; 
Horns et al. 2006). Recently, the discovery of TeV emission 
from G54. 1+0.3 by VERITAS (Acciari et al. 2010) presents 
a brand new case for highlighting the relative significance of 
hadrons in PWNe. 

G54. 1+0.3 is a Crab- like (Lu et al. 2002) SNR with 
properties very similar to the Crab Nebula in both mor- 
phology and photon spectral indices. The central pulsar, 
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PSR J1930+1852, has a period of P = 137 ms and a pe- 
riod derivative of P — 7.5 x 10~ 13 s s _1 , corresponding to 
a current spin-down luminosity of L s d = 1.2 x 10 37 erg s _1 
and a characteristic age r c ~ 2900yr (Camilo et al. 2002). 
A faint X-ray shell was most recently detected surrounding 
the PWN up to ~ 6' from the pulsar (Bocchino, Bandiera, 
& Gelfand 2009). The SNR has been suggested to be at a 
distance of 6.2 kpc by the HI absorption and morphological 
association with a molecular cloud (Leahy et al. 2008). 

Recent AKARI observation discovered an infrared (IR) 
loop, which is explained to be a star-formation loop around 
the G54.1+0.3 PWN (Koo et al. 2008) and is alterna- 
tively explained to be the freshly- formed dust in the su- 
pernova ejecta (Temim et al. 2009). Using VLA radio polar- 
ization and Spitzer mid-IR observations, Lang et al. (2009) 
found a molecular cloud located at the southern edge of 
the PWN and suggested an interaction between the PWN 
and the cloud. In 7-rays, VERITAS observed the VHE TeV 
emission from G54.1+0.3 and found that the efficiency of 
converting the spin-down energy to 7-ray emission is high 
and the ratio of 7-ray to X-ray luminosity is as large as 
0.7. This ratio, two orders of magnitudes higher than that 
of the Crab, is the highest among all the nebulae sup- 
posedly driven by young rotation-powered pulsars (Acciari 
et al. 2010). This may imply that the VHE TeV emis- 
sion has extra components in addition to the contribution 
from commonly-acknowledged energetic leptons scattering 
background photons. The newly discovered IR loop and/or 
molecular cloud around the G54. 1+0.3 PWN may act as 
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an appropriate target for the energetic protons to account 
for high-efficiency 7-ray production from this unusual source 
(Bartko & Bednarek 2008). 

In this letter, we show that the TeV emission from 
G54.1+0.3 cannot be accounted for by leptons alone, but 
can be naturally explained by introduction of a hadronic 
component. 



Thus the normalization parameter Qo(t) can be derived as 
(Tanaka & Takahara 2010) 
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2 MODEL AND RESULTS 
2.1 The Pure-Lepton Case 

We first try to reproduce the wide-range radiation spectrum 
of G54.1+0.3 from radio to TeV using a pure lepton compo- 
nent. For calculating the spectral evolution of the PWN, we 
specify the evolution of the time-dependent injection spec- 
trum and that of the magnetic field in the following. 

Let us consider the relativistic wind of leptons produced 
within the light cylinder of the pulsar where the spin-down 
power L(t) is injected into PWN. A termination shock is 
formed in the outflowing relativistic wind, where the ram 
pressure is balanced by the pressure of surrounding medium, 
and accelerates particles to high energies. The leptons pro- 
duced inside the light cylinder of the pulsar account for the 
radio emission, while the wind leptons accelerated by the 
shock have a Fermi-type energy spectrum and contribute to 
the X-ray emission. 

As usual, we assume that the injection spectrum of the 
relativistic particles Qmj(7,i) obeys a broken power-law 
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where Qo is normalization coefficient, 7 is the Lorentz factor 
of the relativistic electrons and positrons, and the minimum 
(7min), maximum (7 max ), and break (7b) Lorentz factors to- 
gether with the energy indices [p\ and P2) are assumed time- 
independent. Parameter 7 max is obtained so as to confine the 
accelerated electrons within the PWN (i.e., the electrons's 
Larmor radius must be less than the radius of the PWN) 
(Venter & de Jager 2006) 
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where magnetization parameter a is the ratio of the electro- 
magnetic energy flux to the lepton energy flux at the wind 
shock of the PWN. Parameter 7mm = 100 is assumed so as 
to reproduce the flux of the observed minimum frequency at 
radio wavelengths. Bucciantini et al. (2010) found that 7b is 
at a similar value in a narrow range of 10 5 -10 6 for several 
PWNe of a variety of ages, which is closely related to the 
working of pulsar magnetospheres, pair multiplicity, and the 
particle acceleration mechanisms. Therefore, here we adopt 
7b = 5 x 10 5 without loss of generality. 

The injection spectrum can be related to the spin-down 
power L(t) of the pulsar at given time t by assuming that a 
fraction (rj e ) of the spin-down power is converted into lep- 
ton luminosity: rj e L{t) — j Q(j, t)^m e (? dj . For a spin-down 
pulsar, L(t) = L [l + (t/T-o)]" (n+1)/(n " 1) , where L is the ini- 
tial spin-down power, to the characteristic timescale, and n 
the breaking index (here we adopt n = 3 for simplicity). 



On the assumption of magnetic-field energy conserva- 
tion (see Tanaka & Takahara 2010 for the comparison of 
various approximations of magnetic field evolution), 
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the time-varying field strength of the nebula is given by 
1/2 
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where t/b is the fraction of spin-down energy converted to 
the magnetic energy and J?pwn the average radius of the 
PWN. (In parenthesis, the magnetization parameter is thus 
essentially a ~ t]b/t) c .) Because the young G54.1+0.3 PWN 
(~ 2900yr) may be in an evolution stage before the re- 
verse shock passage (typically at 1 x 10 4 yr, e.g., Reynolds 
& Chevalier; Gelfand 2009), we also assume that the PWN 
is freely expanding at velocity upwn and thus have wpwn ~ 
550(7?PWN/1.8pc)(t/2900yr)- 1 km s -1 . 

The volume-integrated particle number as a function of 
energy is described by the continuity equation in the energy 
space: 



J^V( 7 ,t) + — [7( 7l t)JV(7,*)] = 



where 7(7, t) is the cooling rates of the relativistic leptons 
including the synchrotron radiation, the inverse Compton 
scattering off the cosmic microwave background (CMB) and 
ambient IR radiation, and the adiabatic expansion, i.e., 

7(7, *) = 7syn(7, t) + 710(7) + 7ad(7, *), (7) 

and Tesc is the escape timescale and can be estimated as in 
Bohm diffusion (e.g., Zhang et al. 2008), 
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where the current magnetic field strength 80^iG (see j)3.1|l is 
used. The adiabatic loss 7 a( j = — j/t is the dominant cooling 
process for the low energy particles and insignificant for the 
high energy ones. 

The time-dependent lepton distribution is numeri- 
cally solved from the continuity equation (|6}. Then multi- 
wavelength non-thermal emission can be calculated for the 
process of synchrotron radiation and inverse Compton scat- 
tering, with photon spectra plotted in Figures [T] and [2] (as 
described below). 

Here the 7-rays are considered to purely come from lep- 
tons scattering the soft radiation field (CMB, IR and opti- 
cal photons in the Galactic plane, and the IR-optical-UV 
emission of possible young stellar objects (YSOs) in the IR 
loop). The IR background at the Galactic disc is character- 
ized by temperature 25K and energy density two times larger 
than the CMB, while the optical background by tempera- 
tures between 5000 and 10 4 K and energy densities equal 
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Figure 1. Comparison of the predicted spectra in the pure-lepton 
Models A (solid line) and B (dashed line) with the observed data 
for G54.1+0.3 in radio (Natasha et al. 2008; Lang et al. 2009), X- 
rays (Lang et al. 2009) and 7-rays (Acciari et al. 2010). The model 
parameters are described in the text of £|2.1I The red dashed line 
shows the 1 year, 5<r sensitivity for the Fermi LAT (Fermi LAT 
2007). 

to the CMB. The incident IR photons from the SNR are 
defined by a ~ 90 K blackbody radiation with the energy 
density ~ 5.3 x 10~ 12 erg cm" 3 based on the Spitzer IRAC 
fluxes at 24/jm and 70(im from Temim et al. (2009), a fac- 
tor of roughly 5 larger than the IR energy density in the 
Crab Nebula and 13 larger than the energy density in the 
CMB. In the calculation we also take into account the pos- 
sible IR-optical-UV starlight from 11 possible YSOs, which 
has an energy density ~ 4.4 x 10 _11 erg cm" 3 with a black- 
body temperature T ~ 35000-ft: (Koo et al. 2008). The IC 
flux is dominated by scattering with the IR photons from 
the SNR, while the IC scattering with other components 
are insignificant by comparison. Note that the power of 
synchrotron self-Compton emission to synchrotron emission 
Pssc/Psy* = U S yn/U B < 10" 2 (here eq.(27) in Tanaka & 
Takahara 2010 is used), the contribution of 7-ray emission 
for G54. 1+0.3 PWN from IC scattering off the synchrotron 
radiation is negligible. 

For the physical parameters of the PWN, we set Ln « 
1.4 x 10 39 ergs s _1 , 75 = 5 x 10 5 , and pi = 1.2 according to 
previous studies (Camilo et al. 2002; Lang et al. 2009; Buc- 
ciantini et al. 2010) and leave other three parameters, rj e ,riB, 
and P2, adjustable. For comparison, we develop three sets of 
parameters for leptonic model. In Model A, we reproduce the 
observed results of radio to X-ray emission (which are syn- 
chrotron) and get r/ e = 6%, t\b = 8%, and P2 = 2.4. As can 
be seen in Figure [T] (the solid line), the resulting TeV emis- 
sion from leptons is lower than the observed flux by more 
than an order of magnitude. In Model B, we change parame- 
ters to reproduce the observed TeV emission by IC scattering 
soft photon fields described above. The adopted parameters 
are rje = 92%, t\b = 8%, and P2 = 2.1. The resulting syn- 
chrotron radio and X-ray emission (the dashed line in Fig- 
ure [TJ excess the observation data by more than an order of 
magnitude. The current magnetic field strength for Model A 
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Figure 2. The same as Figure [T] but for pure-lepton Model C 
(solid line) . The parameters are described in the text of i|2,ll The 
IC flux (the solid line on the right side) is dominated by scatter- 
ing with the IR photons from the SNR, while the IC scattering 
with the IR photons from Galactic diffusion (dashed), the CMB 
(dotted), and the starlight of the possible YSOs (dashed-dotted) 
are also shown. 

and B, 80/iG (derived from observation, see i )3.ip . has been 
used in Eq.©. In order to match both the synchrotron and 
IC emission to the observed data, we explore the parame- 
ter space and obtain the third model (Model C) (Figure [5} 
with = 99.8%, r]B = 0.15%, and P2 = 2.8. However, this 
corresponds to a weak magnetic field ~ 10/jG. If we adopt 
an age of 2000 yr for this PWN as obtained by Bocchino 
et al.(2009) in their dynamic evolution model, other than 
2900 yr, then lower field strength would be needed in Model 
C. Such low values of the field strength are inconsistent with 
that derived from observation, as will be discussed in £ ]3.1I 
Therefore, it is hard for a pure-lepton model to reproduce 
the radio, X-ray, and TeV data simultaneously, and thus the 
leptons alone cannot account for the 7-ray emission. 

2.2 The Lepton-Hadron Hybrid Case 

We now consider the contribution to the TeV emission from 
a hadronic component besides the leptonic contribution. In 
this model, both leptons and ions extracted from the charged 
polar cap region are accelerated in the rotating magneto- 
spheres of neutron stars and PWN termination shocks (e.g., 
Zhang et al. 2009). For simplicity, we assume the protons 
gain energy from central pulsar and are represented by a 
power-law spectrum which is common for the acceleration 
process. Then the total energy of protons of the PWN is 
W p = fA p E~ ap E p dE p = J Q rj p L(t')dt' , where r\ p is the 
energy fraction converted to protons, A p the normalization 
coefficient and a p the spectral index of accelerated protons. 
So the energy released from the pulsar consists of the kinetic 
energy of particles (rj e and r] p ) and the magnetic energy (t/b)- 
For the energy, E p , of the accelerated protons, the rest en- 
ergy of protons (9.4 x 10 s eV) is adopted as minimum and 
the energy at the "knee" (3x10 eV) as the maximum. Note 
that the energy converted into leptons and magnetic field in 
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Model A is only a small fraction (r/ e = 5% and rjB = 8%, 
respectively) of the total spin-down energy of central pulsar. 
In fact, in the study of the Vela X PWN, Horns et al. (2006) 
have questioned where the remaining energy injected from 
pulsar is and suggested a hadronic origin of TeV emission. 
Hence, we assume r\ v — 87% in the lepton-hadron hybrid 
case (denoted as Model D). 

The Bohm diffusion timescale of the PWN particles de- 
termined from Eq. JSJ (~ 10 4 yr) is much longer than the 
PWN age. Therefore, the protons are considered to be well 
confined in the PWN and the escape losses of protons are 
negligible. The cooling time of p-p interaction is (e.g., Aha- 
ronian 2004) t pp « 1.8 x 10 6 (nb/30cm" 3 ) -1 yr, much longer 
than the age of G54. 1+0.3, where rib is the average density 
of target baryons in the PWN (see below). Hence the colli- 
sion losses of the PWN protons are negligible as well. Also 
because t pp is almost energy-independent in the energy re- 
gion above lGeV, the total spectrum of protons remains un- 
changed (Aharonian 2004). The contribution from the sec- 
ondary leptons that are created by protons interaction to 
the overall spectrum is negligible too, as compared with the 
dominant contribution of the primary leptons (Horns et al. 
2006; Zhang et al. 2009). 

In addition to the contribution from the leptons as given 
in Model A, we calculate that from p-p interaction so as 
to match the observed TeV flux. For the tt° decay ensu- 
ing from p-p collision, the analytic emissivity developed by 
Kelner et al. (2006) is used. It is difficult to determine the 
detail process of energetic protons captured by the bary- 
onic targets, since this process depends on geometry of the 
PWN and the targets and anisotropy of the magnetic field 
and diffusion coefficient. Thus, we assume that a small frac- 
tion (£) of all hadrons is captured by baryonic targets (as 
suggested by Bartko & Bednarek 2008). The wide-range 
spectrum of the PWN can now be well reproduced with 
£ ~ 8 x 10 -3 (nb/30cm -3 ) -1 and the results are shown in 
Figure [3] Here a target baryon density ~ 30 cm -3 has been 
adopted from the estimate of the IR clump density (Temin 
et al. 2010); this number can also be typical of the density of 
the molecular materials, which Koo et al. (2009) and Lang 
et al. (2010) reported to detect. Apparently, even such a low 
capture efficiency is sufficient for hadrons to produce the 
observed flux of TeV emission. 



3 DISCUSSION 

In the pure-lepton case, Model C seems to marginally match 
the wide-range spectrum of the G54.1+0.3 PWN; by com- 
parison, however, the lepton-hadron hybrid case (Model D) 
can reproduce the spectrum better and more physical in the 
following aspects. 

3.1 Magnetic field 

In i]2.1l the field strength obtained in Model C (the lepton 
case) is 10/j.G or even lower. Such values of field strength are 
actually weaker than that derived from observation. Based 
on radio luminosity, Lang et al. (2009) derived an equipar- 
tition field of 38/iG. However, they suggested stronger field 
in the light of the strong polarization which is organized 
on large scales of the nebula and implies the PWN is filled 
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Figure 3. The same as Fig[T] but for lepton-hadron hybrid Model 
D. The parameters are described in the test of £|2.2I The solid 
line on the right side is dominated by 7r° decay ensuing from 
p-p interaction. The inverse Compton scattering with IR photos 
from SNR (dashed-dotted line), IR photos from Galactic diffusion 
(dashed), starlight of the possible YSOs (dashed-dottcd-dotted) 
and the CMB (dotted) are also shown. 



with magnetically-dominated plasma. They also found an 
alternative field strength of 80-200/iG by using the lifetime 
of the X-ray emitting particles. In Model D (the lepton- 
hadron case), however, we use 80^iG which can typify the 
field strength estimated by Lang et al. 



3.2 TeV index 

In Model C, the calculated TeV slope (~ 2.6-3) of the IC 
spectrum cannot well match the VERITAS data point (with 
photon index 2.4, Acciari et al. 2010). Matching the TeV 
slope would entail a lepton ensemble with a unreasonable 
large energy index 3.8. Even if the energy losses in high en- 
ergy leptons are considered, we, using the time-dependent 
model, find the energy index of accelerated leptons by rela- 
tivistic shock is 2.8, still considerably higher than the univer- 
sal power-law index 2.2-2.3 for Fermi-type acceleration by 
the shock of large Lorentz factor using different approaches 
(e.g., Horns et al. 2007). As a contrast, the observed slope 
is easily reproduced by protons p-p interaction with a mild 
proton index a p = 2.4. This proton index is fortuitously 
similar to the lepton index that is used to reproduce the 
synchrotron X-rays in Model A. 



3.3 Baryonic targets 

The IR loop closely around the G54. 1+0.3 PWN discovered 
by AKARI was suggested to be star-forming region (Koo et 
al. 2008), while it was also argued to be the freshly formed 
supernova dust heated by early-type stars belonging to a 
cluster in which the supernova exploded (Temim et al.2009). 
It was also reported that a molecular cloud is found to be 
located at the southern edge of the PWN by the VLA radio 
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and Spitzer mid-IR observations and thus an interaction be- 
tween the PWN with the cloud was suggested (Lang et al. 
2009). These components within or surrounding the PWN, 
whatever they are, may readily be a baryonic target for the 
bombardment of the PWN protons, and therefore it is very 
reasonable to expect the 7-ray contribution from the hadron 
interaction. This is the very case that we address in Model 
D. This scenario seems to naturally explain the exceptionally 
high 7-ray to X-ray luminosity ratio of G54. 1+0.3 among all 
the rotation-powered PWNe. 

The Fermi observation at GeV band will be important 
to discriminate between the leptonic model and the hadronic 
model. In the pure-lepton model (cases A, B, and C; see 
Figures Q] and [2}, the theoretical GeV 7-ray flux of the PWN 
is basically below the 1 year, 5a sensitivity of the Fermi LAT, 
while the lepton-hadron hybrid model (case D; see Figure [3]) 
predicts a GeV flux above the sensitivity. 



4 CONCLUSION 

We have calculated the multi-band non-thermal emission 
from the G54.1+0.3 PWN in both the pure-lepton case and 
the lepton-hadron hybrid case. In the lepton case, we find 
that the leptons that are responsible for the radio and X-ray 
synchrotron cannot alone account for the TeV 7-ray emission 
by IC scattering. An addition of hadron contribution by p- 
p interaction can well reproduce the observation spectrum. 
The lepton-hadron hybrid scenario is strongly supported by 
the most recently discovered IR loop and molecular cloud 
in or closely around the PWN. This scenario can also shed 
light on the study of the PWNe with high 7-ray to X-ray 
luminosity ratios. 
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